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^ i Abstract 

PJ ' We propose a differential difference equation in TZ^ x Z^ and study it by Hirota's bilinear 

r^ ■ method. This equation has a singular continuum limit into a system which admits the reduction 

i_^ , to the Davey-Stewartson equation. The solutions of this discrete DS system are characterized by 

^^ , Casorati and Grammian determinants. Based on the bilinear form of this discrete DS system, 

^ ' we construct the bilinear Backlund transformation which enables us to obtain its Lax pair. 

^: 
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^. ■ 1 Introduction 

'k>( I The nonlinear Schrodinger equation (NLS) 

g^ + (vv)g-^2^ ^^^^^^^^^^ -^^iMg = o (1.1) 

a, 13 



dkadkp dxadxp 



is the simplest universal model for the slow evolution of the envelope g(r, t) of an almost monochro- 
matic wavetrain exp(ikor — iLo{kQ)t) in a weakly nonlinear medium of nonlinear dispersion relation 

Lo{k) = uJo{k) + uJi{k)\q\^ + ■ ■ ■ (1.2) 

In a d + 1 dimensional space this equation has d + 1 canonical forms. For d = 1 they are the 
"self-focusing NLS" 

qt + i{Qxx + \q\'^q) = 0, if ujqUJi > (1.3) 
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and the " self-defocusing NLS" 

qt + i{qxx-\q\'^q) = 0, if wq^i < 0. (1.4) 

In the hterature, one can find many results on eqs. 1)1.31 11.4(1 and their discrete versions (see, e.g. 
^ E] and references therein) . 

There are several important generalizations of the NLS's (|1.3l 11.41) . The best known example 
is the Davey-Stewartson (DS) equation [21ll[ , a partial differential equation in TZ^ given by: 

iqt + cr^iqxx + (^^qyy) = -a\q\'^q + #, (1-5) 

4>xx- (^'^(pyy = '^a{\q\'^)xx, (1-6) 

where a^ = 1,0"^ = ±1. In the hyperbolic case a^ = 1, the system ()1.51 II. 6|) is called DSI while 
DSII corresponds to the elliptic case a^ = —1. In the following we will focus on the DSI equation. 
The DSI equation ()1.5| ll.6|l is a reduction of the system 

iqt + iqxx + qyy) = -aq'^r + qcj), (1.7) 

-in + (r-xx + ryy) = -aqr"^ + r(j), (1.8) 

(pxx - 4>yy = ^a{qr)xx, (1-9) 

obtained by requiring that r = q* . By the variable transformation dx = —m{dx + dY),dy = 
-75 (9x — dy), 4> = cnqr + aij), the system (|1.7m.9() may be transformed into 

iqt + {qxx + qYY) = aq^, (1.10) 

-irt + {rxx + ryy) = ar^, (l-H) 

i^XY = \{dl+dl.){qr). (1.12) 

Integrable discrete versions for the DS equations have not been much studied yet although much 
work has been done on the discrete NLS equation jlj. A few partial results have been presented by 
NijhofF and Konopelchenko in jTOJI^El- The purpose of this paper is to propose a new discrete 
integrable system of equations which can be considered as a discrete version for the DSI system 
(|1.1()H1.11JI) . 

Let us consider the following system 

ivt + aie""-i+""+i-2"r;„_i + a2e'''=-i+"'»+i-2"r;fc+i - (ai + a2)v = 0, (1.13) 

-iwt + aie""-i+''"+i-2«ii,^_^^ + a2e"'=-i+"'=+i"^"u;fc_i - (ai + a2)w = 0, (1.14) 

Zi - 2;ie""+l.'=+l+"-"'=+l-""+l + Z2Vk+lWn+l = 0, (1.15) 

where ai,a2-,zi and Z2 are constants. In eqs. ()1.13tfn3|) and in the following we always use a 
simplified notation for f(n,k,t). We write explicitly a discrete independent variable only when it 
is shifted from its position. For example, 

f = f{n,k,t), fn+i = f{n + l,k,t), fk+i = f{n,k + l,t), f„.+i,k-i = f{n + l,k - l,t). 

Let us now show that eqs. (|1.13tffTT3|l may be thought of as a discrete version of the DSI system 
(|1.1()H1.12B . Let us set 

2 2 1 1 

ai = -;t, a2 = -pr, zi = -p, Z2 = --rO, en = X, 5k = Y, 
e^ 0-^ eb 4 



and expand the dependent fields in power series around 5 = and e = 0, 

Vk+i = v{n, k) = q{ne, {k + 1)5) = q{X, Y + 5) = q + 5qY + -W^YY "^ 

e2 
Wn+i = w{n + \,k) = r{{n + l)e, k5) = r{X + e,Y) = r + erx + -w^xx + 

Un+i = u + eux + -^'^xx H 

Un+i,k+i = u + eux + 5uy + -^uxx + eduxY + -^uyy H 



Then the continuous hmit of eqs. (|1.13ffTTT3|) gives 

iqr + {qxx + qyy) = -2q{d\ + d^)u, (1.16) 

-irr + {rxx + ryv) = -2r(ai + d^)u, (1.17) 

uxY = -jaqr, (1.18) 

where idr = idt — |9x + fcV- Under the transformation V' = —2a{d\ + dY)u, T — *• t, the system 
(|l.l6Hl.l8l) reduces to the DSI system (I1.1()H1.1:JI) . 

In the following we will study eqs. (ll.l.StffTT^ using Hirota bilinear method. By the dependent 
variable transformation 

u = lnF, w = e-*("i+"2)*G/F, w = e'^'^'+'^^^^H/F (1.19) 

eqs. ()1. 13111. l3|) are transformed into the bilinear form 

[iDt + aie-'^" +a2e^'=]G-F = 0, (1.20) 

[iDt + aie-'^" +a2e^'^]F-H = 0, (1.21) 

where, as usual, the bihnear operators Df and exp{6Dn)\9^ are defined as: 

exp(5L'„.) a • b = a{n + 5)b{n — 5). 



Dl'a ■ b 



t'=t 



In section 2, we present the double-Casorati determinant solutions to the differential-difference 
system (|1.13ll.l5p . Its Grammian determinant solutions are presented in section 3, while in section 
4, a bilinear Backlund transformation and Lax pair are derived. Section 5 is devoted to the 
conclusions and a discussion of the result obtained. 



2 Double-Casorati determinant solutions to the discrete DS sys- 
tem 

It is well-known that the continuous DS equation (|1.7H1.9|) has solutions expressed in terms of 
double- Wronskians of the solutions of the Spectral Problem El El El IHI- In this section, we 
present the double-Casorati determinant solutions for the discrete DS system H1.20MTT^^ . An 



example of double- Casorati determinant solution for eqs (|1.13ffTTT^ is the one-soliton solution 
given later in Fig. 1. 

Let us introduce the following double-Casorati determinant: 



|0,1,--- ,m-l; 0',1', 



(2iV-m-l)'| 



02 (ra) 



(f>i{n + m - 1); V'i(fe) 
(/>2(n + m - l)■,^lJ2ik) 



ijJi{k + 2N -m- 1) 
'ilj2lk + 2N - m - 1) 



<t)2N{n) ■■■ 4>2N{n + m-l)]'il)2N{k) ■■■ 'ip2N{k + 2N - m - 1) 
where (j)rin, t) and ^rik, t)(r = 1, 2, • • • , 2A^) satisfy the equations 

d 

i—ipr{n) = -anpr{n - 1), 

d 

^w:ipr{k) = a2'ipr{k- 1). 
ot 

Taking into account eq. (|2.1() we can state the following Proposition: 

Proposition 1 The following double-Casorati determinants 

F = |0,1,--- ,m-l; 0',1',--- , (2iV - m - 1)'|, 
G = zi|0, !,••• ,m; O', !',••• , (2iV - m - 2)'|, 
H = — 10,1,--- ,m-2; 0',l',--- ,(2Af-m)'|, 

Z2 



provide solutions to eqs. M.2U\^1~^) . 

Proof: From eqs. ()2.4H2.6|) for any integer number i and j we have 



n+i,m+j 



|1,2,--- ,m-l + i; 0', !',••• ,{2N - m - I + j)'\, 
Gn+i,k+j = 2:i|0, 1,--- ,m + i; l',2',--- , (2iV - m - 2 + j)'|, 

H,,+i^k+j = -\l,2,--- ,m-2 + i; O', !',••• , (2iV - m + j)'|. 

2^2 



(2.1) 



(2.2) 
(2.3) 



(2.4) 
(2.5) 
(2.6) 



(2.7) 
(2.8) 



Prom the equations (|2.2[ 12. 3j) we get 

iFt = -ail -1,1,--- ,m-l; 0',l',--- ,(27V-m-l)'| 

+a2|0,l,--- ,"1-1; (-1)',1',--- ,(27V-m-l)'|, 

iGt = zi(-ai| -1,1,--- ,m; 0',l',--- ,(27V-m-2)'| 

+a2|0,l,--- ,rn- (-l)',l',--- , (27V - m - 2)'|), 

ii^j = —(-ail -1,1,--- ,m-2; O', l', • • • , (2iV - m)'| 

Z2 

+a2|0,l,--- ,m-2; (-l)',l',--- ,(2iV-m)'|). 

Introducing eqs. H2.7H2.12l) into eq. (11.20(1 we get the determinant identity [HJ: 

|1,2,--- ,m; 0',1',--- ,(2iV-m-l)'||0, 1,--- ,m-l; l',2',--- , (27V - m)'| 

-|1,2,--- ,m; l',2',--- , (2iV - m)'||0, 1, • • • ,m-l; 0',l',--- ,(27V-m-l)'| 
+ |0, 1,--- ,m; l',2',--- ,(2iV-m-l)'||l,2,--- ,m-l; 0',l',--- ,(27V-m)'| 



(2.10) 
(2.11) 

(2.12) 



= 0. 

(2.13) 



Substituting eqs. (|2.7H2.11j) into eqs. (|1.21l[0^ we get the equations 

Qi(| -1,0,--- ,m-2; 0',1',--- , (2iV - m - 1)'||1, 2, • • • ,m-l; O', !',••• ,(27V-m)'| 
-I -1,1,--- ,m-l; 0',1',--- ,(2iV-m-l)'||0, 1,--- ,m-2; 0',l',--- ,(27V-m)'| 
+ 1 -1,1,--- ,m-2; 0',1',--- ,(2iV-m)'||0, 1,--- ,m-l; 0',l',--- , (27V - m - 1)'|) 
+a2(|0,l,--- ,m-l; l',2',--- , (2iV - m)'||0, 1, ■ • • ,m-2; (-1)',0',--- ,(2iV-m-l)'| 
-|0, 1,--- ,m-2; (-1)',1',--- ,(2iV-m)'||0, 1,--- ,m-l; 0',l',--- ,(2iV-m-l)'| 

+ I0,1,--- ,m-l; (-l)',!',--- ,(2iV-m-iy||0,l,--- ,m-2; 0',l',--- , (2iV - m)'|) = 0, 

aid -1,0,- •• ,m-l; 0',1',--- ,(2iV-m-2)'||l,2,--- ,m; O', !',••• ,(2iV-m-l)'| 
-I -1,1,--- ,m; 0',1',--- ,(2iV-m-2)'||0,l,--- ,m-l; 0', !',••• ,(2iV-m-l)'| 
+ |-1,1,--- ,m-l; 0',1',--- ,(2iV-m-l)'||0,l,--- ,m; 0',l',--- , (2iV - m - 2)'|) 
+a2(|0,l,--- ,m; l',2',--- , (27V - m - 1)'||0, 1, • • • ,m-l; (-1)',0',--- ,(2iV-m-2)'| 
-|0, 1,--- ,m-l; (-1)',1',--- ,(2iV-m-l)'||0, 1,--- ,m; 0',l',--- ,(2iV-m-2)'| 

+ |0,1,--- ,m; (-1)',1',--- ,(2iV-m-2)'||0,l,--- ,m-l; 0',l',--- , (2iV - m - 1)'|) = 0, 

which are identically satisfied when we take into account the determinant identities H2.13() . In this 
way Proposition 1 is proved. 

To construct the soliton solution, we choose a simple solution of eqs. (|2.2|l2.;-i|) 

2N 2N 

Mn, t)=Y, ariPT'^e^'^'P'', MK t) = Y. fewC'e-*"^"'*, r = 1, 2, • • • , 2iV (2.14) 

1=1 1=1 

where pr, q^o-r, br are arbitrary constants. Than the one-dromion solution of the discrete DS system 
is obtained by setting N=l in eq. 1)2. 1() and choosing (j)r{n),iJr{k) (r = 1,2) given by eq. ()2.14|) 
with, for example, ai = i, 02 = —i, z\ = Z2 = 1. In such a case we have 

F = {aub2i - aaifoiOpr^r'e"^^^""'^^* + («12&22 - a226i2)P2>2"'e-(P^+''^)* 

+(011622 - a2i6i2)pr>2~'e-(Pi+«^)* + (012621 - a226ii)p2>r'e-(P2+^i)* (2.15) 

G = (aiia22 - a2iai2)pr>2 "(P2"' " Pi')e-^'"-^P'^' (2.16) 

H = (611622 - h2ib^2)qi\2\<i2^ - gr')e-(^i+^^)*. (2.17) 

In Fig. 1, we plot the 1-dromion solution \v\ = ■^, \w\ = ^-pr m. the nk-plane with an = 022 = 
i, ai2 = 0, a2i = 1, 611 = |, 612 = 5, 621 = -\, 622 = 0, pi = e, p2 = e"\ qi = e^, 52 = e"^, at 
the time t = 1. 

3 Grammian determinant solutions to the discrete DS equation 

The Grammian technique was first used by Nakamura for constructing the solutions expressed in 
terms of the special functions for the two-dimensional Toda lattice equation and the KP equation 
|13l I14j . In [7j we can find a Grammian determinant solution for the continuous DS system. In 
this section, we present solutions of the discrete DS system written down in terms of Grammian 
determinants. At the end we show that by a proper choice of parameters the double-Casorati 
determinant solution and Grammian determinant solution give the same 1-soliton solution. 
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Figure 1: The l-dromion solution: (a) l^il — field , {b) {wl — field 



Proposition 2 The functions 



F = \c + n\ = \F\. 



(3.1) 



G = ^ 

«! 



$(n+ 1) 



^'{-k + lf 



if 



2:2 "2 



F ^(/c + l) 

$'(-n+l)^ 



(3.2) 



where F is a {M + N) x (M + N) matrix, C = {c^u) is a (M + A^) x (M + N) constant matrix, Vt 
is a {M + N) X [M + N) block diagonal matrix 



n 



I -00 "PrinW A -n)dt 



f;°° Mkm-k)dt I ' 



and <I>, $', ^, ^' are M + N column vectors 

$(n) = (c^i(n),--- ,ipM{n); 0, • • • ,0)^, 
$'(-n) = ((^;(-n),--- ,V7M(-n); 0, ■ • • ,0)^, 
*(fc) = (0,--- ,0; Mk),--- ,i^N{k)f, 
^'(_A;) = (0,...,0; V'i(-fc),---,V'^(-fc))^, 

wii/i ipr{n,t), if'j{n, t),ips{k,t),'ijj'i{k,t), r,j G {1, • • • ,M}, s, / € {1, ■ ■ ■ , A^}, satisfying the following 
equations: 

dip'j{-n) 



dipr (n) 

' dt ~ 

solve the equations lil.2(An~^} . 



-anpr{n- 1), 
a2ipsik-l), i 



dt 
dt 



-- ai(fj{-n- 1), 
-a2^p'i{-k- 1), 



(3.3) 
(3.4) 



Proof: Using eqs. (|3.3[ 13. 4|) . we are able to express, after some calculations, the functions 
appearing in eqs. ()1.2Ufnr^^ in terms of the Grammian determinants 



F, 



n+l 



F $(n + l) 



ai 



<^'{-ny 



Fk+i = F 



02 



^'{-k^ 



^(A; + l) 




(3.5) 



'k-1 



F + 



Oi2 



F ^(k) 

^'(-k + \Y 



Fn~X = F + 



ai 



$'(-n+l) 



$(n) 




(3.6) 



Fn 



+l,fe+l 



02 



F ^{k + 1) 



^'{-k) 



T 










1 



$(n+ 1) 




aia2 



F $(n+l) *(A: + 1) 

$'(-n)^ 

*'(-A:)^ 



(3.7) 



G 



k+l 






F $(n + l) 

*'(-fe)^ 



^n 



+1 



2^2 "2 



F ^(A; + l) 

<^'{-nf 



(3. 



H, 



fc-i 



^;2a2 



F ^(/c) 

$'(-n + l)'^ 



Gn-l 






F $(n) 

^'(-/c + l)^ 



(3.9) 



F ^{k) 

^'{-kf 



^(n) 



$'( 



-n 



(3.10) 



iG. 



ai 



{-02 



^'(-A:) 



T 



$(n+ 1) 




ai 



^'(-A: + l) 



$(n) 







+ i 



^'(-A: + l)^ 
$'(-n)^ 

F $(n+l) ^(/c) 

^'(-A: + l)^ }, 

"^'{-kf 



<^(n+l) <I>(n) 




(3.11) 



iHt 



Z2a2 



-{«! 



F ^(/fc + 1) 

$'(-n)^ 



02 



F ^(/c) 

$'(-n + l)^ 



^'(-/fc)^ 

F ^(/fc + l) $(n) 

$'(-n + l)^ 

$'(-n)^ 



F ^(/c + 1) ^(A;) 

$'(-n + l)^ 





}■ 



(3.12) 



We can thus prove that the functions F, G and H given by eqs. (|3.11 13.21) effectively satisfy 
the discrete DS system as, by substituting eqs. (|3.5H3.6|) into eqs. (|1.2UffO^ we get the foUowing 
three Jacobi identities for the determinants 



F $(n + l) 

<^'{-nf 



F ^(/fc + 1) 

^'{-kf 

F $(n + l) 

^'{-kf 



F 



F $(n + l) ^(fc + 1) 

$'(-n)^ 

^'(-fc)^ 







0, 



(3.13) 



F $(n) 

$'(-n+l)^ 



F ^(/c + l) 

$'(-n)^ 



+ F 



F $(n) 

<^'{-nf 



F 



F ^(A;+l) $(n) 

$'(-n + l)^ 

<P'{-n)^ 







} 



F ^(/t) 

^'{-kf 



F ^(A; + l) 

$'(-n + l)'^ 



F ^(A; + l) ^(/c) 

$'(-n + l)'^ 

'^'{-kf 



'^'(-A;)^ 



^(A;+l) 




$'(-n+l 



\T 







} = 0. 



(3.14) 



^'(-A;)^ 



$(n+ 1) 



^'(-A; + l) 



T 



^(A:) 







+ F 



F $(n + l) ^(A;) 

^'(-A; + l)^ 





F $(n+l) 

*'(-A: + lF 



F ^(A;) 



} 



+ { 



F $(n + l) 

^'{-k+lf 



F $(n) 

$'(-n)'^ 



F $(n+l) 

$'(-n)^ 



F $(n + l) $(n 

1''(-A; + 1)'^ 

$'(-n)^ 

<l>(n 



^'(-A: + l 



\T 







} = 0. 



(3.15) 



The simplest soliton solution for the discrete DS system ()1.13ffr!T3|) is obtained by taking the 
simplest possible choice for the functions ipr-,^'j-,'4^s-,'4^'i satisfying eqs. ()3.31l3.4p . i.e. an exponential 



n Aa\kr t 



^r{n) = Ke'""' 



k „—ia20J, t 



where ki,kj,u;s,u}i are arbitrary constants. 
When N = 1, if we take 






c 



-^ 

Ql 



ki = uJi = 2, ki = {- + i) -^, uji = (- + i) ^, 



we have the following 1-soliton solution for equations (|1.13() - ()1.15|) : 

u = ln(^— [(1 + 2i)"(l + 2i)'=e("i-"2)* - 
aia2 

fc-l it[-iai + (i-i)a2] 



1]), 



V = a2Zi 



2^+Hl + i) 



w 



ai 



[1 + 2i)"(l + 2i)^e(°i-"2)* - 1 

2^+1(1 + ^\n-l ^it[{^-i)ai + ^a2] 



Z2 (1 + 2i)'^(l + 2i)'=e("i-"2)* - 1 ' 



(3.16) 
(3.17) 
(3.18) 



This same solution is obtained by considering the double-Casorati determinant solution ()2.4H2.6|) 
with A^ = 1 and 

4 Bilinear Backlund transformation and Lax pair 

In this section we construct a bilinear Backlund transformation for the bilinear equations (|1.2fll - 
ll.22() . and then we derive from it a Lax pair for the discrete DS system (|1.13MTTT^ . 

To do so, let us redefine the functions F, G and H in term of one function / depending on an 
additional discrete variable m 

F(n, k; t) = f{m, n, k; t), G{n, k; t) = f{m + 1, n, k] t), H(n, k; t) = f{m — 1, n, k] t), 

Then eqs. (fTMfT^ can be written as: 



[zie 



l/2{D„~Dk) 



+ z^e'mD,-D^)+D,^ ^ ^^gi/2(D„+D,)]^ . ^ ^ o_ (4 2) 



We can now state the following proposition: 

Proposition 3 The bilinear system i4-l[ \4-^ ^^^ ^^c Backlund transformation 

[^^gl/2(D^+D,) _ g-l/2(D,„+D,) _ ^^^II2(D^-D,)^f . ^ = Q, (4.4) 

[iDt - ai^e^'"-^" - as^e-^'^]/ -5 = 0, (4.5) 

Pi P2 

where Pi, P2, fJ-i, M2 o'^e arbitrary constants, with fii, ^2 satisfying the constraint 

mzi + ^2^2 = 0. (4.6) 



Proof: Let / be a solution of equations (|4.HI4.2j) . Using eqs. (|4.;-iH4.6|) . we can by straightfor- 



ward calculations show that eqs. (|4.11 14.2|) are satisfied for g{m, n, k; t) 

= 2sinh(l/2i?™)(iA/ • 5) • /5 + 2ai sinh(l/2(A - DrnW^^""^ f ■ g) ■ (e-^/^^'"/ • g) 
+ 2a2 sinh(l/2A)(ei/2(Dfe+D,„) J . ^^ . (g-i/2(D,+D„)_^ . ^ 

= 2sinh(l/2D„)(iA/ • 5) • /9 + 2ai sinh(l/2(A - Z)„))(^e^™-i/2^"/ • g) • (e^^/^D^/ • 5) 

Pi 

+ 2a, sinh(l/2A)(^eV2(^™-^.)/ . 5) . (e-i/2(Z).+D.)^ . ^ 

P2 

= 2sinh(l/2D„)(iA/ ■g)-fg- 2ai^ sinh(l/2Z)^))(e^'"-^"/ • 5) • /5 

Pi 

- 2a2^ sinh(l/2D„))(e-^'=/ • g) • /g = 0, 

P2 

= 2zi sinh(-l/2A)(eV2^"/ • 5) • (e-V^^-/ • 5) 

+ 2z2sinh(l/2(An - A))(e^/'(^-+^^-)/-5) • (e-i/2(D„+^'=)/ • 5) 
= -2Z1/.1 sinh(l/2A)(ei/2^"/ • g) ■ (e^'"-i/2D„^ . ^) 

+ 2Z2M2 sinh(l/2(Z?„ - A)(ei/2(^™+^'=)/ • 5) • (ei/2{i?^-i?.)/ . g) 
= -2zi^i sinh(l/2A)(ei/2^"/ • g) • (e^'"-i/2i?n^ . ^) 

- 2Z2M2 sinh(l/2A)(ei/2^"/ • g) ■ (e^™-i/20./ . ^) = q 

In this way, Proposition 3 is satisfied and eqs. (|4.3H4.6j) constitute a BT for (|4.11I4.2|) . 

From the bilinear Backlund transformation (|4.3l4.6j) . we can derive a Lax pair for the discrete 
DS system (|l.l3lll.lKI) . 
Let us set 

u = \nf, v = ——, w = ^—, 4>=-;- (4.7) 

Under the dependent variable transformation ()4.7p . the bilinear BT ()4.3l4.5p become the nonlinear 
equations: 

M - 4>n+l - ^^lVWn+l4>m-l,n+l = 0, (4.8) 

l32Wk-l(t)m-l,k-l - Wk-l(t> - tJ-2W(j)m-l = 0, (4.9) 

icPt + ai^«,_iu;„+ie""+i+""--2>„_i.„+i + a2^e"'=+i+"'=-i-2>,+i = 0, (4.10) 

Pi P2 

where /3i, f],, /"i, /j., are arbitrary constants satisfying the constraint (|4.6I) . Eliminating 0m-i,n+i, 4'm-i,k-i, <t>m-i 
from eqs. (|4.8H4.l(7]l . we obtain the following Lax pair for the differential-difference DS system 

r. (Pl<t)n-l,k-l-(t>k-l. /3l(/)„-l -(/> 

P2( ■ )-^2( ) -(t)Wk-l=Q, (4.11) 

y"l1'n-l,fc-l fJ-lVn-1 

Pi l^ P2 
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By imposing the compatibility of eqs. 1)4.111 I4.12|l we obtain the discrete Davey-Stewartson system 
(|1.13l4TTT3|) . In fact, from eq. H4.11() . we can derive 

Pl4>n-l,k-l = (f>k-l + -^Vn-l^k-lifpWk-l -\ {Mn-l " 0)], (4.13) 

P2 fJ'lVn-1 

(f)n+l,k-l = Mk-l - -r'^k~l[Wn+l,k~l(pn+l H {M " (f^n+l)], (4.14) 

P2 ^J^■lV 

II R 

Mn-l,k+l = ^k+1 -\ Vn-l,k+l[ (/3l0n-l " 0) " (j^k+iw], (4-15) 

the expres sions of 4)n~i,k-i, (pn+i,k-i, (pn-i,k+i in terms of (/>, (fi n-i, (pk-i, (p n+i, (pk+i- By differenti- 
ating eq. ()4.11|) with respect to t and substituting into it eqs. H4.13tH?T3|l we obtain an expression 
in terms of just (p,(pn-i, 4'k-ii (pn+i and (pk+i- Equating to zero the coefficients of cp , (/>n-i, 4>k-ii 
(pn+i, and (pk+i we derive that the coefficient of (pn-i gives eq. (|1.13p . the coefficient of cj) gives eq. 
()1.14|) . both the coefficients of (pn+i and (pk+i give eq. ()1.15|) and the coefficient of (pk-i vanishes. 

5 Conclusion. 

A discrete version of the Davey-Stewartson (DSI) system is proposed and investigated using the 
bihnear method. This DSI system exhibits N-sohton solutions expressed in terms of determinants of 
two different types, double-Casorati and Grammians. Moreover, we have constructed the bilinear 
Backlund transformation and derived from it its Lax pair. 

A few problems are still open. Among them the most significant is surely to find the proper 
reduction which gives the Davey-Stewartson equation from the system. Moreover, since in the 
continuous case we have two physically interesting cases, the DSI and DSII equations, it would also 
be interesting to find the discrete version of the DSII equation. 
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